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0\ • Abstract 

• Most boson emitting sources contain a core of finite dimensions surrounded 

,1^ ■ by a large halo, due to long-lived resonances like io,r],r]' , etc. When the 

•'-H . Bose-Einstein correlation (BEC) function of the core can be determined we 

X 

. show that its intercept (A) measures, as a function of momentum, the square 

of the fraction of core particles produced. A simultaneos measurement of BEC 
and the single-particle distributions can thus determine the characteristics of 
the core. 

There are two types of scales present in the space-like as well as in the time- 
like directions in a model-class describing a cylindrically symmetric, finite, 
three-dimensionally expanding boson source. One type of the scales is related 
to the finite lifetime and geometrical size of the system, the other type is 
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goverened by the change of the local momentum distribution in temporal and 
spatial direction. If the geometrical sizes of the core are sufficiently large the 
parameters of the BEC function obey the mt-scaling observed in SPb and 
PbPb reactions at CERN. The model can describe the measurements of the 
single- and two-particle distributions in the central region of SPb reactions. A 
fit to experimental data shows that the freeze-out of hadrons occurs at a larger 
volume and at a much lower temperature than that given by the measurement 
of the inverse slope of the mrspectrum and standard BEC analysis. 

I. INTRODUCTION. 

Hanbury-Twiss correlations were discovered 40 years ago, revealing information about 
the angular diameters of distant stars The method, also referred to as intensity interfer- 
ometry, extracts information from the quantum statistical correlation function for (partially) 
chaotic fields. The method has extensively been applied to the study of the freeze-out ge- 
ometry in high energy nucleus-nucleus collisions as well as in elementary particle reactions. 
For recent reviews see refs. |@J^. 

Accumulating evidence indicates that the space-time structure of the pion emission in 
heavy ion reactions at the 200 AGeV bombarding energy region at CERN SPS has a peculiar 
feature, namely that the boson emission can be approximately divided into two parts: the 
centre and the halo The central part corresponds to a direct production mechanism e.g. 

hydrodynamical evolution or particle production from excited strings, followed by subsequent 
rescattering of the particles, while the surrounding halo corresponds to pions emitted from 
the decay of long-lived hadronic resonances, like uj,t], rj' and K^, which have a mean decay 
length of more than 20 fm. These long-lived resonances will decay and give rise to mainly low 
momentum pions. However, the momentum distribution of the emitted pions is independent 
from the position of the decay which in turn is smeared out over a large region due to the 
long mean decay length. This mechanism will result in pions of similar (low) momentum 
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from a very large volume. This is in contrast to the core which when expanding in a 
hydrodynamical manner will have pions of similar momentum from a relatively small volume. 
We will investigate the phenomenological consequences of such a structure where there is 
a characteristic length-scale of the core which is smaller than any length-scale that can be 
connected to the surrounding halo. In the case of halo length-scales larger than 20 fm they 
give rise to a sharp peak of the correlation function in the Q < h/R ^ 10 MeV/c region. This 
region is strongly influenced by the Coulomb interaction so the accuracy of the correlation 
function determination will be dependent on proper Coulomb corrections which constitute a 
problem in itself. Furthermore, this Q < 10 MeV/c region is very difficult to measure due to 
that we must properly determine two very close tracks which make experimental systematic 
errors largest in this region of the correlation function even for high resolution experiments 
like the NA35 and NA44 high energy heavy ion experiments at CERN. 

We will show that even if we cannot measure reliably the correlation function below 
Q < 10 MeV/c important new information about the particle production in a core/halo 
scenario can be given by analyzing both the invariant momentum distribution of single 
particles and the correlation function for Q > 10 MeV/c. We will apply our approach to 
NA44 data on SPb reactions at 200 AGeV/c. 

We investigate the core/halo scenario utilizing analytical results in the Wigner-function 
formalism. This formalism is well suited to describe both the invariant momentum distri- 
bution (IMD) and the Bose- Einstein correlation function (BECF). We utilize the version 
of the formalism discussed in refs. and applied to analytic calculations in refs. [p!0| -[T5 
recently. 



II. THE WIGNER-FUNCTION FORMALISM 

In the Wigner-function formalism, the one-boson emission is characterized by the 
emission function, S{x,p), which is the time derivative of the single-particle Wigner 
function. Here x = (t,r) = {t,rx,ry,rz) denotes the four-vector in space-time, and 
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p — {E,p) — {E,px,Py,Pz) stands for the four-momentum of the on-shell particles with 



mass m — — p^. The transverse mass is defined as mt = Jm? +p1+Py and 



y = 0.5 log 



^ E + p^ 



stands for the rapidity. Attention is drawn to the point that the 

\E-Pz) 

normal size character x stands for a four- vector, normalsize y is the symbol for the rapidity 
while the subscripts ^ and y index directions in coordinate and momentum space. 

For chaotic sources the two-particle Wigner-function can be expressed in terms of the 
symmetrized products of the single-particle Wigner functions. Based on this property, the 
two-particle Bose-Einstein correlation functions can be determined from the single-particle 
emission function solely. These relations are especially simple in terms of the Fourier- 
transformed emission function, 

S{M, K)^ j d^x S{x, K) exp(iAA; • x), (1) 



where 



Pi +P2 

Ak^pi-P2, (2) 



and Ak ■ x stands for the inner product of the four-vectors. The invariant momentum 
distribution of the emitted particles is given by 

d'^n 

^S{Ak^Q,K ^p). (3) 



dydm\ 

The IMD is normalized to unity as 



r d^n 

/ dydm] = 1. (4) 

J 7 7 9 



dy dml 



The two-particle BECF-s are prescribed in terms of our auxiliary function, 

I S{Ak,K) |2 

C{K,Ak)^l + , (5) 



as was presented e.g. in ref. @J^. The effect of final state Coulomb and Yukawa interactions 
shall be neglected as implicitly assumed by the above equation, and we assume completely 
chaotic emission. 

The Wigner-functions are in general complex valued functions, being the quantum ana- 
logue of the classical phase-space distribution functions. For the purpose of Monte-Carlo 
simulations of the Bose-Einstein correlation functions the off-shell Wigner-functions are ap- 
proximated by the on-shell classical emission functions from the simulation in case of 
analytic calculations they are usually modeled by some positive valued functions |To|-[l5|. 
Neither of these approximations seems to be necessary: in fact we need only simplifying 
assumptions as hsted below. 



Assumption 0. There is a resolvable length-scale in the core 

Assumption 1. The bosons are emitted either from a central core or from the surrounding 
halo. The respective emission functions are indicated by Sc{x,p) and Sh{x,p), and fc indi- 
cates the fraction of the bosons emitted from the central part. According to this assumption, 
the complete emission function can be written as 



Both the emission function of the core and that of the halo are normalized similarly to the 
complete emission function, normalized by eqs. (|l|J^,^). 

Assumption 2. We assume that the emission function which characterizes the halo 
changes on a scale Rh which is larger than Rmax ~ f^/Qmin, the maximum length-scale 
resolvable by the intensity interferometry microscope. Qmin is taken as 10 MeV/c in our 
considerations. However, the smaller central part of size Rc is assumed to be resolvable. 



III. THE CORE/HALO MODEL 



S{x,p) = fcSc{x,p) + (1 - fc)Sh{x,p). 



(6) 



Rff > i?< 



max 



> R, 



(7) 
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This inequality is assumed to be satisfied by all characteristic scales in the halo and in the 



central part, e.g. in case the side, out or longitudinal components |T6| , p!0| of the correlation 
function are not identical. 

Based on Assumption 1 the IMD can be re- written as 

= fc + (l-/c) , (8) 
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dy dm^ dy dm^ dy dm^ 

where the subscripts c, h index the contributions by the central and the halo parts to the 
IMD given respectively by 

d'^Ui 

= S,iAk = 0,K = p) (9) 



dy dmi 

for i = c,h. The normalization condition according to Assumption 1 reads as 

d^rii 

dydm] = 1 (10) 

dy drn^ 

for i = c,h. Note that the IMD as expressed by eq. (H) includes the possibility that the 
IMD for the bosons of the halo is different from the bosons emerging from the central part. 
Thus the relative contribution of the halo and the core is a function of the momentum in 
this model. 

The BECF is expressed by 

\fJ,{Ak,K) + {l-f,)S^X^k,K) |2 
C{K, Ak) = l + , (11) 

I fJ,iAk = 0,K = p) + il-f,)ShiAk = 0,K = p) p 

which includes interference terms for boson pairs of (c, c) (c, h) and {h, h) type. Due to the 
assumption that the emission is completely chaotic, the exact value of the BECF at the 
Ak = is always 2 in this model. 

The measured two-particle BECF is determined for Q > Qmin ~ 10 MeV/c, and any 
structure within the Ak < Qmin region cannot be resolved. However, the (c, h) and {h, h) 
type boson pairs create a narrow peak in the BECF exactly in this Ak region according to 
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eq. (|TT]), which cannot be resolved due to Assumption 2. According to this assumption, the 
Fourier-transformed emission function of the halo for non-zero relative momenta vanishes at 
the given resolution Qmm- At zero relative momentum the same quantity gives the single 
particle momentum distribution from the halo, which is not effected by the two-particle 
resolution. 

Including the finite resolution effect, symbolized by the horizontal bar, the measured 
BECF can be written as 

\Sc{Ak,K)\^ 

C{Ak, K) = l + K{K=p- Qrmn) , (12) 

I Sc{Ak = 0,K = p)\'^ 

where 

I fcSc{Ak = 0,K = p)\'^ 

X^{K = p]Q,^in) = , (13) 

I /, S.iAk = 0,K = p) + il- /,) ShiAk = 0,K = p)\' 

or expressed with the IMDs 

\fc P 

dy dmf 

K{p;Qmin) = • (14) 

I |2 



dy dm^ 

The IMD is connected to the cross-sections by the relations 

1 (i^cTj d'^Ui 

=< rii > (15) 

a dy dml dy dm^ 

< nc > 

for i = c,h. From Eq. (8) we deduce fc = . The effective intercept K{p;Qmin) is 

<n> 

then also given by 

d'^o'c 

I |2 



dy dml 



K{p]Qmin) = • (16) 



d^ 



a 

|2 



dy dmi 



The effective intercept parameter A* shall in general depend on the mean momentum 
of the observed boson pair, and we have < A* < 1. The equation above shows that this 
effective A* has a very simple interpretation as the square of the fraction of core particles 
to all particles emitted. We also see that this equation makes it possible to measure the 
cross-section of the bosons produced in the central core with the help of the combined use 
of the BECF parameter A* {y, nit) and the cross-section of all particles. 

Note that the A*(?/, m^) intercept parameter in our picture reflects the drop in the BECF 
from its exact value of 2 due to the contribution from the resonance halo, when measured with 
the finite two-particle momentum resolution Qmin- Thus the measured intercept, X^{y,mt) 
does not coincide with the exact value of the BECF at zero relative momentum within this 
description. 

The effect of the halo is to introduce a momentum-dependent effective intercept param- 
eter to the BECF. The shape of the BECF for Q > Qmin shall be solely determined by the 
freeze-out phase-space distribution in the central part. This central part or core is usually 
well accessible to hydrodynamical calculations. 

Note also that this result has in principle nothing to do with the similar relation obtained 
for the case of partially coherent, partially chaotic fields, where A* = fl^^, and fine =< 
nine > / < n-tot > gives the fraction of the total multiplicity in the chaotic field, although 
the relation formally is similar. For the case of partially coherent fields, the Bose-Einstein 
correlation function contains an interference term in between the coherent and the chaotic 
fields which leads to a double-Gaussian or double-exponential structure [jl^, emphasized 



recently in ref. |]T8[. In our case, the BECF contains not three but only two terms, and it 
looks very much like the BECF would look like for the case the halo were missing i.e. fc = 1. 

Note also, that the focus is not on the reduction of the intercept of the BECF due to the 
halo (resonance decays) but more precisely on the momentum- dependence of this reduction. 
The fact of the reduction has been known before however its momentum depen- 

dence was not utilized as far as we know. Within our formalism however it is straightforward 
to show that the core fraction of bosons can be estimated by 

8 



/oo fOO I d n 
dy / dmlJX{y,mt) (17) 
-oo Jin^ , , 9 



dy dm^ 
and we also have 



min 



Kv^^t) < fc < meLX^JX{y,mt). (18) 



Application. Present NA44 data []T9[ for central S + Pb reactions at CERN SPS with 
200 AGeV bombarding energy show an approximately independent intercept parameter: 
= 0.56 ± 0.02 and 0.55 ± 0.02 at the quite different mean transverse momenta of 150 
MeV and 450 MeV, respectively. This suggests that in the midrapidity region the momentum 
distribution of the resonance halo is similar to that of the central part, and the halo contains 
1 - = 25 ± 2 % of all the pions. For this special case of A = = const we have the 
following simplified equations for the IMD and the BECF: 

d'^n d'^Uc d^Uh 

(19) 



dy dml dy dm^ dy dm^ 



C{Ak,K) = l + f^R,{Ak,K). (20) 

I.e. the only apparent effect of the halo is to reduce the intercept parameter of the measured 
BECF to A = /2 while the IMD and the (Afc, K) dependence of the BECF is determined 
by the central part exclusively. 

Pions coming from the resonance-halo surrounding a strongly interacting centre are pre- 
dominant at lower values of the transverse momentum according to the SPACER calculations 
in ref. [0] and RQMD calculations in ref. 0, similarly to the results of calculations with 



HYLANDER and the hydrodynamical calculations of the Regensburg group [^. The 
resonance fraction as a function of the transverse momentum was explicitly shown in the 
publications 1^,0. Each of the mentioned models predicts that the halo of resonances pro- 
duces predominantly low momentum pions, and they predict that the effective intercept 
K{p',Qmin) should increase with increasing transverse mass in contrast to the measured 
values. 
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The constancy of the effective parameter A=k (y, rrit) with respect to rrit suggests a mech- 
anism of enhancing the low momentum particles in the core. Such an enhancement has a 
natural explanation in a hydrodynamical description of an expanding core as we will present 
below, see also [^,|T^. The mechanism is a simple volume effect, as pointed out in [0. To 
calculate the cross-section you integrate over the volume from where the particles of a given 
momentum are emitted. This volume is measured by the BECF. In the NA44 experiment 
all three radii of this volume are equal and they are inversely proportional the the square 
root of rrit HH. This will give an extra volume factor, oc (m^)"'^/^, enhancing the single 



particle cross-section at low rrit. We will return to this point. 

Discussion. The general result for the BECF of systems with large halo coincides with 
the most frequently applied phenomenological parameterizations of the BECF in high energy 
heavy ion as well as in high energy particle reactions 0. Previously, this form has received 
a lot of criticism from the theoretical side, claiming that it is in disagreement with quantum 



statistics |]T8[ or that the A parameter is just a kind of fudge parameter or a measure of our 
ignorance, which has been introduced to make theoretical predictions comparable to data. 
Now one can see that this type of parameterizations can be derived with a standard inclusion 
of quantum statistical effects - if we assume that we discuss interferometry for systems with 
large halo. The fact that most of the studied reactions, including + e~ annihilations, 
various lepton-hadron and hadron-hadron reactions, nucleon-nucleus and nucleus - nucleus 
collisions are phenomenologically well describable [0] with variations of eq. (12) thus does 
not exclude the possibility that these high energy reactions all create boson emitting systems 
which include a large halo. 

The applicability of the halo picture to a given reaction is not necessarily the only one 
possible explanation of a reduced intercept. It is known that the final state interactions 



may have an influence on the effective intercept of the two-particle correlation functions ||22 
and the effect has been shown to result in stronger drop in the intercept value for smaller 
source i.e. for hadronic strings created in lepton-lepton, lepton-hadron or hadron-hadron 
collisions. However, more detailed calculations indicate that the higher order corrections for 
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the final state interactions may very well cancel the first order effects resulting in a very 
small total final state interaction correction [^]. In our calculations these effects have been 
neglected, so our results should in principle be compared to data which present the genuine 
Bose-Einstein correlation function. This function is not necessarily the same as the short- 
range part of the two-particle correlation function since Coulomb and Yukawa interactions 
as well as the rj'ri decay chain and other short-range correlations may mask the quantum 
statistical correlation functions. However, as the particle density is increased (in the limit 
of very large energy or large colliding nuclei), the Bose-Einstein correlations together with 
the final state interactions dominate over other short-range correlations due to combinatoric 
reasons. 



IV. THE MODEL OF THE CORE 



For central heavy ion collisions at high energies the beam or z axis becomes a symmetry 
axis. Since the initial state of the reaction is axially symmetric and the equations of motion 
do not break this pattern, the final state must be axially symmetric too. However, in order 
to generate the thermal length-scales in the transverse directions, the flow-field must be 
either three-dimensional, or the temperature must have significant gradients in the transverse 
directions. Furthermore, the local temperature may change during the the duration of the 
particle emission either because of the re-heating of the system caused by the hadronization 
and/or intensive rescattering processes or the local temperature may decrease because of 
the expansion and the emission of the most energetic particles from the interaction region. 

We model the emission function of the core for high energy heavy ion reactions as as 



Sc{x, K) d'^x 
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I 



I 



vni cosh(r7 — y) exp 



K ■ u{x\ 



^i{x] 



uix) 



cosh(?7) 



V 



1 + b' 



r + r' 



2r, 



/ 



b — , b — , smh.{ri) 

To To 



H{T)dT Todrj dr^ dvy, (21) 



»,2 _j_ 2 ^ ^ 

,2 y 



2t, 



/ 



(22) 
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T 

n^) = ^4 r (23) 

1 + 



^ 7-2 _|_ 2 ^ 

r, X ' ' y 

1 + 

V J 



V J 



= . (24) 

T(x) To 2Rl 2A7f 



We include a finite duration, H{t) oc exp(— (r — ro)^/(2Ar^)), r = yW^^z^, tq is the mean 
emission time, Ar is the duration of the emission in (proper) time, u{x) is the four-velocity 
of the expanding matter, /i(x) is the chemical potential and T{x) is the local temperature 
characterizing the particle emission. 

This emission function corresponds to a Boltzmann approximation to the local momen- 
tum distribution of a longitudinally expanding finite system which expands into the trans- 
verse directions with a non-relativistic transverse flow. The decrease of the temperature in 
the transverse direction is controlled by the parameter a, while the strength of the transverse 
flow is controlled by the parameter h. The parameter c = 1 is reserved to denote the speed 
of light, and the parameter d controls the strength of the change of the local temperature 
during the course of particle emission. 

For the case ofa = 6 = (i = 0we recover the case of longitudinally expanding finite 
systems. The finite geometrical and temporal length-scales are represented by the transverse 
geometrical size Rg-, the geometrical width of the space-time rapidity distribution Ar] and by 
the mean duration of the particle emission Ar. We assume here that the finite geometrical 
and temporal scales as well as the transverse radius and proper-time dependence of the 
inverse of the local temperature can be represented by the mean and the variance of the 
respective variables i.e. we apply a Gaussian approximation, corresponding to the forms 
listed above, in order to get analytically trackable results. We have first proposed the 
a = 0, 6 = 1 and = version of the present model, and elaborated also the a = h = d = Q 
model |T2] corresponding to longitudinally expanding finite systems with a costant freeze- 



out temperature and no transverse flow. Soon the parameter h has been introduced p6[ and 
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it has been realized that the transverse flow has to be non-relativistic at the saddle-piont 
corresponding to the maximum of the emission function. Yu. Sinyukov and collaborators 
classified the various classes of the ultra-relativistic transverse flows |^5|], and introduced 
a parameter which controls the transverse temperature profile, corresponding to the a = 
case. We have studied PT| the model-class a, 6, d = which we extend here to the d ^ Q 



case too. 

The integrals of the emission function are evaluated using the saddle-point method 
2^^6|. The saddle-point equations are solved in the LCMS |12|, the longitudinally co- 



moving system, for r/^ << 1 and r^^s « tq- These approximations are warranted if 
\y ~ l/ol << 1 + ^if"f^tlTQ and [3t = Pt/rrit « (a^ + The flow is non-relativistic 

at the saddle-point if Pt << (o^ + The radius parameters are evaluated here up 

to O(r^ s/ro) -|- 0{ris), keeping only the leading order terms in the LCMS. However, when 
evaluating the invariant momentum distribution (IMD), sub- leading terms coming from the 
cosh.{ri — y) pre-factor are also summed up, since this factor influences the IMD in the lower 
rrit region where data are accurate. 

A. Bose-Einstein correlations 

The experimental BECF is parameterized in the form of C {Q l,Q side, Q out) = 1 + 
Aexp(— i^lQl — R^sideQlide ~ R-outQout) where the intercept parameter A and the radius pa- 
rameters may depend on the rapidity and the transverse mass of the pair. The parameters 
of the the correlation function are given by 

R-lide = R-l, (25) 
i?L = R' + f^T^rl (26) 
Rl = r^^vl (27) 



We obtain [|12],|T3| 
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Ill 

— = — + — cosh[77,] (28) 

r>2 d2 d2 
-n-* J^G ^ 

111 1 

= + cosh[77,] , (29) 

/\r]l /\rf ArjT cosh^[7?s] 

111 

= + cosh^ [77,] . (30) 



Ar^ At' At| 
where the thermal length-scales are given by 

i?|=^^ — , (31) 

a' + b' rrit 

Ar]^ = —, (32) 

mt 

Ar| = — — . (33) 
(f rrit 

As a consequance, the parameters of the BECF-s are dominated by the smaller of the 
geometrical and the thermal scales not only in the spatial but in the temporal directions 
too. These analytic expressions indicate that the BECF views only a part of the space-time 
volume of the expanding systems, which implies that even a complete measurement of the 
parameters of the BECF as a function of the mean momentum K may not be sufficient to 
determine uniquely the underlying phase-space distribution. We also can see that for pairs 
with I yo — Y \ « 1 + Arj'^Mt/To the terms arising from the non- vanishing values of rjg can be 
neglected. In this approximation, the cross-term generating hyperbolic mixing angle r/^ ~ 
thus we find the leading order LCMS result: 

C{Q; K) = l + X^K) expi-RlQl - RI,,,QI^ - i^L^L), (34) 

with a vanishing out-long cross-term, Rout,L = 0. K{K) is given by eqs. (14) or (16). 
According to the previous part, X*{K) — in case of the NA44 measurements. 
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The difference of the side and the out radius parameters is dominated by the hfetime- 
parameter At*. Thus vanishing difference in between the -Rgj^^e can be generated 

dynamically in the case when the duration of the particle emission is large, but the thermal 
duration Att becomes sufficiently small. This in turn can be associated with intensive 
changes in the local temperature distribution during the course of the particle emission. If 
the finite source sizes are large compared to the thermal length-scales we have 

1 To 

At^ = At| ^ (35) 

At|. (f rrit 
1 + 

Ar^ 

Rl ^ tI — (36) 
mt 

^Le ^ (37) 

+ 6^ mt 

Thus if (i^ >> + 6^ ~ 1 the model features a dynamically generated vanishing difference 
between the side and out radii. 

If the vanishing duration parameter is generated dynamically, the model predicts an mt 
- scaling for the life-time parameter as 

d2 p2 -1 

Ar^ = ^ — , (38) 

Pt 

Note that this prediction could be checked experimentally if the error bars of the measured 
radius parameters could be decreased to such a level that the difference between the out and 
the side radius parameters could be significant. 

If the finite source sizes are large compared to the thermal length-scales and if we also 
have -|- 6^ ~ 1, one obtains an mt -scaling for the parameters of the BECF, 

To (a^ + b^) 1 

Rlide-Rlut-Rl-ro—, valid for A« ~ -. (39) 

mt 6' 
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Note that this relation is independent of the particle type and has been seen in the recent 



NA44 data [|T^. This mj-scaling may be valid to arbitrarily large transverse masses with 
f3t ~ 1 if << 1- Thus, to generate the vanishing difference between the side and out 
radius and the rrit scaling simultaneously, the parameters have to satisfy 6^ << + 6^ ^ 
1 << (P, i.e. the fastest process is the cooling, the next dominant process within this 
phenomenological picture has to be the development of the transverse temperature profile 
and finally the transverse flow shall be relatively weak. 



B. Invariant momentum distributions 



The IMD plays a complementary role to the measured Bose-Einstein correlation function 
TTHl3[. Namely, the width of the rapidity distribution at a given rrit as well as X,, the 



effective temperature at a mid-rapidity i/q shall be dominated by the longer of the thermal 
and geometrical length-scales. Thus a simultaneous analysis of the Bose-Einstein correlation 
function and the IMD may reveal information both on the temperature and flow profiles and 
on the geometrical sizes. E.g. the following relations hold: 



A?/(mf)^ = Arj"^ + Arj'^^mt), and 



1 



/ 



+ 



1-/ 



(40) 



To + Tcirut = m) Tq 

The geometrical contribution to the effective temperature is given by Tg = TR%/F^ and 
the fraction / is defined as / = 6^/(a^ + 6^), satisfying < / < 1. The saddle-point sits at 
ris = {yo - y)/{l + Ar/VAr^^), r^,^ = ptbRl / {toAt]^) , ry^s = and = tq. 

For the considered model, the invariant momentum distribution of the core can be cal- 
culated as 
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exp 



dydmf (Svr)" 



A^o , At^ 

m, (27rAr7,V2)i/2 (2vri?2) cosh(r7,) exp(+Ar/2/2) x 



X exp 



{y - yof 



At 



\ 2iAri' + Ari'^) J 



\ 


mt 


( 






( 


exp 




1 


-f- 




exp 


/ 


[ To 


I 


2j 


) 


V 



2 \ 



2{T^ + Tg)I 



■ (41) 
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This IMD has a rich structure: it features both a rapidity-independent and a rapidity- 
dependent low-pt enhancement as well as a high-pt enhancement or decrease. 

V. APPLICATIONS TO NA44 DATA 

We have fitted the NA44 preliminary invariant momentum distributions for pions and 
kaons together with the final NA44 data [|l9l for the dependence of the BECF parameters 
for both pions and kaons. As A has been shown not to be dependent on rrit the IMDs 
measured are proportional to the IMDs of the core. Fixing the parameter ?/o = 3 we get a 
description of the IMD and the BECF for both pions and kaons. The result of the preliminary 
analysis, see Table, indicates large geometrical source sizes for the pions, Raij^) ~ 7 fm and 
Rl ~ 20 fm, late freeze-out times of tq ~ 7 fm/c, a vanishing duration Ar for both pions 
and kaons and finally surprisingly low freeze-out temperatures, Tq ~ 100 MeV. The kaons 
appear from a smaller region of Rg{K) ^ 4 fm. 

Let us mention, that the measured points for the IMD contain statistical errors only, 
systematic errors are not included. This implies that the errors of the parameters, as stated 
in the Table, are probably to be enlarged, when we can use final data. The Ax^ is relatively 
large from the pion invariant momentum distribution. We have scanned the parts which 
contribute to the large increase in and concluded that this increase seems to be related 
to some large fluctuations of data points at the edge of the acceptance region. 

VI. CONCLUSIONS 

In summary we have studied the case when the central boson-emitting region is sur- 
rounded by a large halo, which also emits bosons. If the size of the halo is so large that it 
cannot be resolved in Bose-Einstein correlation measurements, lot of information shall be 
concentrated in the momentum dependence of the intercept parameter of the correlation 
function. We have shown that with the help of the Bose-Einstein correlation measurement, 
the invariant momentum distribution can be measured for the two independent components 
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belonging to the core and the halo, respectively. The results do not depend on any particular 
parametrization of the core nor of the halo. 

Analysis of the NA44 data for two-pion correlations indicated that the normalized in- 
variant momentum distribution of the pions from the halo of long-lived resonances within 
errors coincides with the normalized invariant momentum distribution of the pions from 
the central core. This result can be explained with a hydrodynamical model of the core 
development in space-time. The number of pions coming from the halo region, > 20 fm, was 
found to be 25 ± 2 % of the total number of pions within the NA44 acceptance. 

Instead of observing a small fireball we are observing a big and expanding snow-flurry 
at CERN SPS S + Pb reactions accordingly to this analysis of the partly prehminary NA44 
data. The central temperature in the snow-flurry is slightly higher than in the outer regions, 
however the temperature gradient is rather small (a^ = 0.03±0.01). The temporal changes of 
the local temperature during the particle emission seem to be rather intensive {(P — 0.9±0.3). 

Acknowledgments: I want to express my sincere gratitude to the organizers of this Sym- 
posium for creating such a nice and stimulating atmosphere. 
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Table 

Results of a simultaneous fit to NA44 preliminary invariant momentum distributions and 
Bose-Einstein radius parameters for tt and K in central S + Pb reactions at 200 AGeV 



X^/ndf, full fit 


490/220 


Ax^, pion singles (128) 


351 


Ax^, kaon singles(94) 


137 


Ax^, radii measurements(9) 


2 


To [MeV] 


100±2 


To [fm/c] 


6.9±0.2 


-RgCtt) [fm] 


6.9±0.4 


Rg,(k) [fm] 


4.1±0.4 


Rg,l [fm] 


21±6 


Ar(7r) [fm/c] 


O.liO.l 


At{K) [fm/c] 


6±4 




0.03±0.01 


6^ 


0.85±0.01 


(P 


0.9±0.3 
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